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M. Poincare visits Jefferson Lab: Relativistic Models of Few-Nucleon 
Systems 
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I discuss relativistic models of few-nucleon systems, with particular emphasis on calculations 
of electron-deuteron scattering and the comparison of these calculations with recent data from 
Jefferson Lab. 



1. Introduction 

In the "standard model of nuclear physics" nuclear properties are calculated using a 
Schrodinger equation in which the degrees of freedom are nucleons interacting via energy- 
independent potentials. This picture can now be solved essentially exactly for light nu- 
clei jll'2j . However, it appears it only describes non- relativistic systems, since the Poincare 
algebra is satisfied approximately: order-by-order in an expansion in momenta over the 
nucleon mass |3|4j . 

Here I review attempts to build a phenomenology of light nuclei which satisfies the 
Poincare algebra exactly, and so is manifestly applicable when few-nucleon systems are 
probed at GeV-scale momentum transfers — the kinematic domain accessed at Jefferson 
Lab. I will focus on elastic electron-deuteron scattering, since it involves the simplest 
non-trivial nucleus, and yet also requires the calculation of iViV-system wave functions 
away from the two-body centre-of-mass frame. 

In Sectional write down the Poincare algebra, and sketch some of the different means 
by which model-builders have attempted to satisfy it. Then in Section |3] I describe the 
calculation of matrix elements for the interaction of the probe (electron) with a relativistic 
bound state (the deuteron). I enumerate the places in which "relativistic effects" can 
occur in the calculation, and show one explicit example of such an effect. I then compare 
a number of relativistic models of ed scattering, and address the issue of how we should 
interpret their widely- varying predictions. I conclude in Section 0] 



2. What is a relativistic model? 

The laws of quantum mechanics are the same in all inertial frames of reference. In 
non-relativistic quantum mechanics the generators of the transformations relating one 
frame of reference to another are the total momentum P (translations), the total angular 
momentum J (rotations), and the boost operator Kg. This last operator can be defined 
by its action on state vectors: suppose Abigail and Beatrice have their coordinate systems 
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aligned at t = 0, but Abigail moves with velocity v relative to Beatrice: 

(x\m)A = (x + vt^t))* = (x|e iKG - v |^(t)) B . (1) 

Considering the combined effect on a state vector of all possible pairs of translations, 
rotations, and boosts then leads to the following commutation relations: 

[P\ P j ] = 0; [J\ P] = ie ijk J k ; [J\ K G ] = ie ijk K k ; (2) 

[J\ P 1 } = ie ilk P k ; [Kq, K g ] = 0; [K G , P>] = -i5 i] M. (3) 

Meanwhile, the invariance of Schrodinger equation dynamics means that the Hamiltonian, 
H, commutes with J and P. H does not commute with Kg: 

[H,P i ]=0- [H,J i )=0; [H,K i G ]=iP i . (4) 

The relations define the Galilei algebra, which relates observations in different 

frames of reference, and must be obeyed in non-relativistic quantum mechanics. 

The Poincare algebra also relates observations in different frames, but now it implements 
Einsteinian, rather than Galilean, relativity. The changes in the algebra seem simple 
enough: the only commutators which change are the last two in (jBJ). They become: 

[K\ K j ] = -ie ijk J k ; [K\P j ] = -i5 lj H. (5) 

Any "relativistic" quantum theory should obey the Poincare algebra. The first of the two 
relations (jSJ states that relativistic boosts K, in contrast to their Galilean cousins, are 
not commutative. The second is the key commutator which makes relativistic quantum 
theories difficult to construct: it encodes the absence of absolute time in relativistic 
theories, since it connects the generator of time translations, H, to the generators of 
spatial boosts and translations. In any interacting theory H is the sum of a free part H° 
and a potential energy V, so Eq. ©(b) means that either K or P must depend on V. 

This is one way to understand Dirac's famous "three forms" of relativistic quantum 
mechanics 0. In instant form (which will be the form I mainly discuss here) J and P 
depend on H° alone, and so are the same in the interacting theory as they are for free 
particles. In contrast, K depends on V, and so will be different for systems with different 
dynamics. We say that the angular and linear momentum are kinematical generators, 
while boosts are "dynamical" . In point form P is dynamical, and K is kinematic. Finally, 
in front form linear combinations of the 9 Poincare algebra commutators are taken, such 
that only two, rather than three, generators need involve V. In front form rotations about 
the x and y axes are dynamical. It has, though, the great advantage that boosts in the 
z-direction are kinematic. 

One example of theories that manifestly obey the Poincare algebra are relativistic quan- 
tum field theories (QFT), such as Quantum Electrodynamics. Field theories can be quan- 
tized so that they use instant-form, front-form, or point-form dynamics l . If instant form 
is adopted then J and P can be written down directly from the stress-energy-momentum 
tensor of the free theory, and do not connect different Fock-space sectors of the theory 2 . 

1 See Ref. jS] for calculations of electron-deuteron scattering in a front-form QFT-based approach. 
2 1 do not consider problems of gauge invariance in such definitions here. 
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In contrast, K inolves the particle- number-changing interaction piece of the QFT Hamil- 
tonian, and so connects different sectors of the Fock space. In a weak-coupling theory 
perturbative forms of the boost are useful, but in nuclear physics, we deal with coupling 
constants of order 1, and some non-perturbative approximation to the boost, which man- 
ages to preserve — at least approximately — the Poincare algebra, is needed. This means 
that the application of QFT to relativistic nuclear physics is not straightforward. 

Manifest covariance is a sufficient, but not necessary, condition for a theory to be 
relativistic. The Bakamjian-Thomas (BT) construction [7j is a technique by which a P, 
J, and K that obey the Poincare algebra may be constructed, once a potential energy V 
is given. Crucially, particle number is conserved in the BT construction — as long as V 
does not change it. In instant form, we first introduce the intrinsic spin j, the operator 
X = «Vp, and the invariant mass of the system M. Some algebra then shows that 
defining [S|: 

l r „„, PxJ 



H = y/M 2 + P2; J=j + XxP; K = -~{H, X} - — — , (6) 

gives a set of operators that obeys the Poincare algebra. Note that K depends on V, but 
the BT P and J are the same as in the free theory. BT constructions for front form and 
point form are also possible |5|. Some results for ed scattering employing them will be 
displayed below. 

Suppose then that we have a relativistic quantum theory, with dynamics — an H° and a 
V — which is God-given in the rest frame of the many-body quantum system. Considering, 
for the time being, only one state in the spectrum, we write: 

# |V)o = Eo\tf>) . (7) 

Relativistic covariance imposes constraints on matrix elements in this theory. For instance, 
if J is an operator corresponding to an observable which is a Lorentz scalar, then we 
should have: 

oWo|V>>o = qWqWq, (8) 

where the subscripts indicate operators and states constructed by observers in frames in 
which the system is, respectively, at rest, and moving with total momentum Q. Quantum- 
mechanical equivalences of the type (JBJ) imply the presence of unitary transformations 
relating wave functions and operators: 

l^>Q = W>o; jQ = U Q JoU j Q . (9) 

Uq is a unitary representation of the Poincare group. Of course, in a manifestly covariant 
theory (jHJ) is satisfied automatically, as long as |-?/>)q and J^q are constructed consistently. 

3. Electron-deuteron as an example 
3.1. Observables 

The ed differential cross section, can (up to corrections of 0(a 2 )) be written in terms 
of two structure functions, A and B: 



da da 
dVt dVt Mott 



A{Q 2 ) + B{Q 2 ) tan 2 (?± 



(10) 
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where 8 e is the c.m.-frame electron scattering angle, q 2 = {p' e —p e 



-Q 2 is the virtuality 



of the photon exchanged between the electron and the nucleus, and ^ Mott is the Mott ed 
cross-section. A and B are related to the deuteron form factors Gc, Gq, and Gm by: 



A = G 2 C + ~r)G 2 M + °-v 2 M*G 2 Q 



B=~ V (1 + V )G 



2 . 

Mi 



V = 



\ 



AM. 



2 ■ 



(11) 



We also consider the tensor-polarization observable, T20, which is defined from the ratios: 



2 G Q 

X = ^Gc 



y = - 3 v 
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~ + (l+v) tan 2 (I 



(12) 
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l + 2(x 2 + y) ' 



(13) 



Measurements of T 20 , A, and B at a fixed Q 2 enable the extraction of Gc, Gq, and Gm, 
as is done, for instance in Ref. |Hj out to Q 2 ~ 2 GeV 2 . In the Breit frame (see Fig.^) the 
deuteron charge, quadrupole, and magnetic form factors are linear combinations of matrix 
elements of the NN four-current J 11 between deuteron magnetic substates, for example: 



G 



c 











(14) 



q/2 



-q/2 



Figure 1. Three momenta of the deuteron and virtual photon in the Breit frame. 



3.2. Calculations 

Thus to calculate A, B, and T20 we either need an explicit calculation of the wave 
functions for deuteron magnetic substates \m z )±q/2, or we must employ a unitary repre- 
sentation of the Poincare group to write: 

q/2(^|^Breitl^)-q/2 = (V> I ^q/2 J Breit C/ -q/2 l^)o- (15) 

To simplify matters, we now expand J^reit = iBreit + ^2B Breit' where is the single- nucleon 
current, and drop J^. This leads to a "Relativistic Impulse Approximation" (RIA). 

If, as we are doing here, we take the view that the Poincare group relates processes 
in different frames, and remain agnostic about relativity's impact on \ip)o, there are two 
potential sources of relativistic effects in the RIA. First, the one-body current: 



JBrcit _ U fi l) 



F 1 (Q 2 ) 1 ^ + F 2 (Q 2 



2M 



w(ki) (16) 



contains effects traditionally called "relativistic", since they arise in a Dirac equation 
treatment of the nucleon, e.g. Thomas precession, the Foldy contribution to (r 2 ). 
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Second, there are boost effects arising from the C/'s. I will mention one simple example 
of these, that the boost incorporates length contraction. In instant form this has the 
consequence that the deuteron wave function, computed in the NN rest frame, is eval- 
uated at a value of p which is length-contracted with respect to the value in the Breit 
frame jlOlllll^j : 

<p[C/- q/2 |'0>o « 4=<P«Ma = -^>(Pc), (17) 

where p c 's components perpendicular to q are the same as p's, but the component parallel 
to q is p • q/rj. Keeping only this length-contraction effect, at present, and approximating 
u^°u by the first term in its p/M expansion — 1 — it is easy to show that: 

q /2<m,| J°K>-q/2 = F 1 (Q 2 ) J ^y 3 r M (p C + <Mp c ). (18) 

This is exactly the result for this matrix element in the non-relativistic impulse ap- 
proximation (NRIA), but with q replaced by q/77. Expanding 77 in powers of Q/M one 
might think that the difference between the NRIA and the instant-form RIA will be very 
small unless Q 2 is large. In fact though, this effect is non-negligible, especially near the 
minimum of Gq at Q ~ 800 MeV, since the charge form factor varies rapidly there. 

3.3. Survey of relativistic approaches 

This then is one relativistic effect, in one approach to relativistic dynamics. What of 
other approaches? A detailed description of various models which attempt to deal with 
this regime can be found in the recent reviews |13I14|15] . I will provide no more than a 
rough sketch here, in each case discussing only the relativistic impulse approximation for 
each model. My plots for these results are taken from the Gross and Gilman review . 

Considering calculations where quantum mechanical wave functions generated in the 
rest frame are boosted to the Breit frame, we can choose instant form, front form, or point 
form. Examples of all three appear in Fig. |21 For the instant form we have a calculation 
by Forest, Schiavilla, and Riska (FSR) [T2j. It proceeds essentially as described in the 
previous subsection, although no p/M expansion is made for the nucleon current, and 
Wigner rotations are included (approximately) in the boost. The input wave function 
used is the AV18. This result includes some two-body-current contributions. 

For the front form we have two representatives: a calculation by Lev, Pace, and Salme 
(LPS), which works in a frame where Q± = 0. In Ref. several input wave functions 
were used, but results displayed here are for the Nijmll wave function. We also show a 
light-front calculation by Carbonell and Karmanov (CK), in which an attempt to restore 
manifest rotational invariance is made ^7]. This is done by allowing the orientation of 
the light front to be an additional vector in the theory, and then eliminating dependence 
on it. Unfortunately, only approximate calculations in this approach exist as yet. 

Moving to the point form, I first note that for Eq. (|16p to be correct I should have 
depending on k' x — ki. In the point form (and also in some versions of the front form |16j ) 
this need not be the same as q [TBI- A consequence of this is that the squared momentum 
transfer to the struck nucleon is f{p)Q 2 i] 2 , where the function / > 1, but depends on the 
nucleon relative momentum p. A calculation by Allen, Klink, and Polyzou (AKP) using 
this form of dynamics, together with the AV18 wave function is displayed below 
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In any one of these forms of quantum mechanics one can take the one-body current op- 
erator and evaluate its matrix element between deuteron wave functions, defining that as 
the relativistic impulse approximation. However, in the front form with Q + = (and also 
in the point-form calculation of AKP) an impulse approximation such as this implicitly 
includes some two-body current effects. That is because the full consequences of Poincare 
covariance and current conservation are not necessarily respected in calculations of indi- 
vidual matrix elements, and must be imposed on the calculation by a procedure equivalent 
to introducing two-body currents which act specifically to restore these symmetries. 

Also shown in Fig. El are two calculations in which |V ; )±q/2 is calculated dynamically, 
that is to say, the dynamical equation of the theory in the moving frame is solved to obtain 
l^). Obviously in order to do this one needs an NN model in which V can be calculated 
in an arbitrary frame, and this is simplest for meson-exchange models. The calculation of 
van Orden, Devine and Gross [THj (vOG) uses an interaction developed in the "spectator 
formalism" and fit to NN data [20]. That of Phillips, Wallace, and Devine (PWD) 
employs a two-body equation which incorporates the effects of Z-graphs and relativistic 
kinematics in the NN system, but only has approximate boost covariance [21] ■ There 
the Bonn-B interaction is chosen. Both calculations have impulse approximations which 
respect deuteron four-current conservation. 

3.4. Results 

Turning now to the results shown in Fig. |2] we see that there are a wide range of RIA 
predictions. All calculations predict a rapid fall in A, but they agree with the experimental 
data to differing degrees, especially once Q 2 ~ 2 GeV 2 . Agreement with the data at the 
level of the RIA is not essential for success, since explicit two-body currents can, and 
probably should, be added to all these calculations. However, large two-body currents 
might cast doubt upon the efficiency of an expansion in hadronic degrees of freedom. For 
further discrimination between models it is good to divide out the rapid overall trend, as 
is done in Ref. [T3]. I will not do a detailed comparison of models with experiment here. 

What is clear is that even were the data to be removed from the slide, there would be 
no definite prediction for the RIA value of A at, say, Q 2 = 2 GeV 2 . This is true even 
though most of these calculations respect the Poincare algebra and current conservation, 
and remains true even if we restrict ourselves to those calculations using the same rest- 
frame dynamics (e.g. AKP and FSR). Symmetries and rest-frame dynamics alone are not 
enough to unambiguously predict A at these high values of Q 2 . 

This is perhaps even clearer in B, where the position of the minimum is an extremely 
sensitive test of the dynamics: both the implementation of the Poincare algebra and 
the choice of potential. But again, it is worth stressing that even calculations with the 
same rest-frame dynamics produce markedly different predictions for the minimum: the 
position of the minimum in the AKP and FSR calculations differs by about 30%. And 
front form and point form even predict opposite directions of the "relativistic shift" of 
this minimum |16|18|22] . 

Finally, in T 20 , most models seem to do reasonably well. The JLab Hall C exper- 
iment j2H|, which extended measurements of t 20 out to Q 2 = 1.7 GeV 2 , shows that 
relativistic hadronic models of deuterium have some success in predicting this observ- 
able. However, detailed comparison with the accurate data at Q 2 < 0.5 GeV 2 reveals 
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Figure 2. Comparison of different relativistic calculations for A, B, and T20. All cal- 
culations use single- nucleon form factors due to Mergell et al. [23 . The calculations, 
in order of their minima in B, are: CK (long dot-dashed line), PWD (dashed double- 
dotted line), AKP (short dot-dashed line), two versions of the vOG calculation (solid line 
and long-dashed line), LPS (dotted line), a quark-model calculation not discussed here 
(widely-spaced dotted line), FSR (medium-dashed line), and a calculation using a p/M 
expansion by Arenhoevel, Ritz, and Wilbois [21], also not discussed here (short-dashed 
line). Figure from Ref. [Tl] . thanks to Ron Gilman and Franz Gross. 
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that models have some problems in precisely reproducing this data. 

To further understand the range of predictions for A and B above 1 GeV 2 I want to 
consider a gedankencalculation, or actually two gedankencalculations. Calculations 1 and 
2 both begin with the same rest-frame and the same one-body j M , but they have two 
different implementations of the Poincare algebra, i.e. 

IV^^I^o; \j>)$ = ug\1>)o. (19) 

Since the two U's are both unitary representations of the Poincare algebra there is a 
unitary transformation that can make an exact equivalence between matrix elements in 
formulation 1 and those in formulation 2. If I choose: 

7(2) _ f/ (2) f/ (l)t (1) (2) t 

then: 

t)MJfM)%2 = { ?/MjM)%2- (2i) 

Crucially though, the unitary transformation ()20|) results in having two-body pieces, 
even though is only a one-body operator. Thus, the equivalence (j21j) only holds if these 
two-body currents are included in . If the relativistic impulse approximation is used 
in formulation 2 also, then the results of the two different calculations will differ, by an 
amount given by the size of the two-body piece induced in Eq. ()2(J|). 

4. Conclusion 

Henri Poincare once said "Les faits ne parlent pas" : facts do not speak. Nevertheless, 
let me attempt to interpret the results presented here. 

I draw two lessons from them. Firstly, there is no unique relativistic impulse approxima- 
tion prediction for electron- deuteron scattering. The RIA is only defined once a particular 
implementation of the Poincare algebra is chosen. Associatedly, the size of the variation 
in Fig. El is indicative of the size of two-body contributions to A and B. 

Does this mean then, that the situation is hopeless? Not at all: it is just that when 
looking at an RIA result one must know how the Poincare algebra was implemented. It 
also suggests that in any given formalism it is important to examine two-body currents. 
Different formalisms will have different two-body currents, and it would seem that in some 
approaches these currents are quite sizeable. Work in this direction has already been done 
in instant form [12], and is progressing in front form [22], and in a calculation involving 
dynamical wave functions [2*o] . 

More data would be useful in order to help different theorists pin down these contribu- 
tions in their formalisms. A theory-independent statement about the size of two-body cur- 
rents is impossible at Q 2 ~ 2 GeV 2 , but once a specific theory is chosen data is needed to 
calibrate the two-body currents which contribute significantly to elastic electron-deuteron 
scattering at these momentum transfers. Especially useful would be more data around 
the minimum of B, since this is a barometer sensitive to many model elements. 

Having worked hard to tune-up our calculations to reproduce elastic electron-deuteron 
data, an important future step is the extension of these calculations to deuteron elec- 
trodisintegration (for steps in this direction see |27|28j ). Finally, much of what was said 
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here also applies to the three-nucleon system. The boosted NN t-matrix is input to the 
Faddeev equations for NNN bound and scattering states. Different implementations of 
the Poincare algebra lead to different boosted NN t-matrices. The unitary equivalence 
between these implementations is maintained only if the associated three-body forces are 
kept in the three-nucleon calculation. 
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